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Abstract 

The QCD treatment of the parton structure of virtual photons is briefly re- 
called, and possible limitations and open questions are pointed out. Various 
models for these densities are compared, completed by a short discussion of 
the treatment of heavy flavors. Finally, different ways to measure the parton 
distributions of virtual photons in e + e~ and ep experiments are summed up. 



1 Introduction 

Theoretical studies of photonic parton distributions of real, i.e., on-shell, photons have a 
long history initiated by Witten's work [[ij. On the experimental side the past few years 
have seen much progress since the advent of HERA. The observation of 'resolved' photon 
induced ep processes, like (di-)jet photoproduction, allows for tests of the hadronic nature 
of (real) photons which are complementary to structure function measurements in e + e~ 
collisions, where new results from LEP/LEP2 have improved our knowledge as well 0. 

Studies of the transition of the (di-)jet cross section from the photoproduction to 
the deep-inelastic scattering (DIS) region at HERA point to the existence of a parton 
content also for virtual photons |3], || . These measurements have revived the theoretical 
interest in this subject and have triggered a series of analyses of the dependence of the ep 
jet production cross section on the virtuality of the exchanged photon || ||. Recently, a 
next-to-leading order (NLO) QCD calculation of the (di-)jet rate in ep (and ej) scattering, 
which properly includes the contributions due to resolved virtual photons, has become 
available |7], || , and resolved virtual photons have been included for the first time also in 
the Monte Carlo event generator RAPGAP ||. 

Pioneering work on the parton structure of virtual photons has been already performed 
a long time ago [10-12]. However, phenomenological models for these distributions have 
been proposed only in recent years [13-16] in view of the expected experimental progress. 
Ongoing measurements at HERA and future structure function measurements at LEP2 
should seriously challenge these models and hopefully lead to a better understanding of the 
transition between the photoproduction and the DIS regime. To finish this introductory 
prelude, let us stress that photons provide us with a unique opportunity to investigate its 
parton content in a continuous range of masses (virtualities) in contrast to the situation 
with nucleons or pions. 

The framework for parton distributions of virtual photons, theoretical expectations, 
and open questions are briefly recalled in Sec. 2. The various different models for the 



parton content of virtual photons are compared in Sec. 3, supplemented by a short dis- 
cussion of the treatment of heavy flavors in Sec. 4. In Sec. 5 we sum up the different ways 
to measure the parton densities of virtual photons in ep and e + e~ experiments. 

2 Theoretical framework: definitions, expectations, 
open questions 

For clarity we henceforth denote the probed target photon with virtuality P 2 = —p 2 by 
7(P 2 ), where p 1 is the four momentum of the photon emitted from, say, an electron in 
an e + e~ or ep collider^]. For real (P 2 = 0) photons we further simplify the notations by 
setting, as usual, 7 = -y(P 2 = 0). 

The concept of photon structure functions for real and virtual (transverse) photons can 
be defined and understood, in close analogy to deep-inelastic lepton-nucleon scattering, 
via the subprocess 7*(Q 2 )7(P 2 ) — > X, as in e + e~~ — ► e ± X ('single tag') or e + e~ — > e + e~X 
('double tag'). The relevant 'single tag' differential cross section can be expressed as in 
the hadronic case in terms of the common scaling variables x and y 



d 2 a(e-f(P 2 ) -> eX) 2na 2 S, 



dxdy Q 4 



(1 + (1 - y) 2 ) F« p2 \x, Q 2 ) - y 2 F2 {p2 \x, Q 2 ^ 



with ^ denoting the photonic structure functions. The measured e + e~ cross section 
is obtained by convoluting ([[]) with the photon flux for the target photon 7(-P 2 ) |17J • The 
range of photon 'masses' (virtualities) produced is 

m 2 e y 2 /(l -y)< P 2 min < P 2 < P 2 max < | (1 - y)(l - cos9 max ) , (2) 

where y is the energy fraction taken by the photon (y = E y /E e ), S is the available 
squared c.m.s. energy, and Q ma x is the maximal scattering angle of the electron in this 
frame. P^ in max in (0) are further determined by detector specifications and/or an eventual 
tagging of the outgoing electron at the photon producing vertex. P^ in effectively measures 
to a good approximation the dominant photon virtuality involved, just as P^ in = m 2 y/{\ — 
y) ~ represents quasi-real photons. Even in the latter case there is, however, still 
a small contribution from the high-P 2 , virtual photon tail of the spectrum, which has 



to be estimated [I6|, |1^] when one tries to extract the parton densities of real photons. 
For transverse virtual target photons j(P 2 ), whose virtuality P 2 is essentially given by 
P 2 ~ Pmini one expects Jit], |TT) a parton content f^ p \x,Q 2 ) along similar lines as for 
real photons. The range of applicability of this 'picture', however, deserves a further 
scrutiny. 

For real photons 7 it is well-known that in the framework of the quark parton model 
(QPM) the x- and Q 2 - dependence of F% L = F^ ^ is fully calculable from the 'pointlike' 
QED process 7*(Q 2 )7 — > qq if one introduces quark masses m q to regulate the mass sin- 
gularities due to P 2 = However, this description is subject to perturbative QCD 



corrections due to gluon radiation not present in the QPM [[T], |2(|. The logarithmically 

1 In the latter case it is common to use Q 2 = —q 2 instead of P 2 , but we prefer P 2 according to 
the original notation used in e + e~ annihilations [jl0| [Tl| | , where it refers to the virtuality of the probed 
(virtual) target photon, and Q 2 is reserved for the highly virtual probe photon 7*(Q 2 ), Q 2 = —q 2 3> P 2 . 







enhanced contributions a s \nQ 2 / 'Qq can be resummed to all orders, removing the depen- 
dence on effective quark masses, where Qq denotes some a priori not fixed renormalization 
point somewhere in the perturbative region Q 3> Aq C d- Of course, this is not the whole 
story, since the photon can undergo a transition into a vector meson of the same quantum 
numbers, which is afterwards probed by the 7*(Q 2 ) (Vector Meson Dominance (VMD) 
assumption). This non-perturbative part obeys the same evolution equations as known 
from the hadronic case. 

Turning to virtual photons, i.e., P 2 7^ 0, it is expected [|H], [Ll| that for large enough 



virtualities P 2 one ends up with a fully perturbative prediction irrespective of Q 2 . To 
facilitate the discussions, it is useful to define the relevant different ranges of P 2 : 

(I) P 2 «A^ CD «Q 2 , (II) P 2 ~K 2 QCD , 
(III) A 2 3CD «P 2 «g 2 , (IV) P 2 ^Q 2 . 

Case (I) we have already discussed above, since it refers to a (quasi-)real photon with 
P 2 ~ 0. In range (III) one can apply similar considerations as long as one restricts oneself 
to transverse virtual photons [10, 11], with the important distinction that P 2 is now 
within the perturbative domain and hence can serve to fix Q , i.e., Qq = 0(V). This is 
the basis for the above mentioned conjecture of absolute predictability in this case, since 
any non-perturbative VMD-inspired contributions are expected to vanish like (1/P 2 ) 2 due 
to such a 'dipole' suppression factor in the vector meson 'propagator'. 

Several questions have to be addressed: up to which values of P 2 (and x, Q 2 ) is the 
non-perturbative part relevant? What are the lower and upper bounds on P 2 in (III), 
i.e., where and how takes the transition to regions (II) and (IV), respectively, place, and 
down to which value of P 2 in (III) should one trust perturbation theory? For smaller P 2 , 
i.e., for a transition to the parton content of real photons (I), one has to find some appro- 
priate, physically motivated prescription which smoothly extrapolates through region (II), 
where perturbation theory cannot be applied, down to P 2 = 0. On the other side, P 2 is 
bounded from above by P 2 <C Q 2 in order to avoid power-like (possibly higher twist) terms 
(P 2 /Q 2 ) n which should spoil the dominance of the resummed logarithmic contributions 
~ a s lnQ 2 / P 2 and, furthermore, to guarantee the dominance of the transverse photon 
contributions in physical cross sections. For P 2 approaching Q 2 (region (IV)) the e + e~ 
result should reduce to the one given by the full fixed order box 7*(Q 2 )t(P 2 ) — > qq in- 
cluding all (P 2 /Q 2 ) n terms and possibly 0(a s ) QCD corrections, which are unfortunately 
unknown so far. 

The question of when fixed order perturbation theory becomes the more reliable pre- 
scription and the concept of virtual transverse photonic parton distributions (i.e., re- 
summations) becomes irrelevant and perhaps misleading is in some sense similar to the 
question of whether heavy quarks should be treated as massless partons or not, which was 
extensively discussed in the literature recently ^2|. Both issues are characterized by 



the appearance of at least two different, large scales, P 2 and Q 2 (or m 2 and Q 2 ), which 
might be indicative for resummations or not. In our case here, however, one is also in- 
terested in the transition to a region where resummations are indispensable (i.e., for real 
photons), but the range of applicability of this approach with respect to P 2 (and possibly 
x and Q 2 ) cannot be determined reliably so far unless the full NLO corrections to the 
l*(Q 2 )l(P 2 ) box will be available to analyze its perturbative stability. 

As already mentioned, for a given Q 2 ^> P 2 and increasing P 2 one expects that the 
resummed results approach the QPM result determined for m 2 <C P 2 <C Q 2 , due to 



the shrinkage of the evolution length, i.e., less gluon radiation. The QPM result can be 
obtained from the process 7*(Q 2 )7(-P 2 ) — > qq, but now P 2 ^ can act as the regulator 
and no quark masses have to be introduced. Taking the limit P 2 /Q 2 — > whenever 
possible, one obtains for F 2 7(p2) [0, |TT| 

^ia^ 2 )=^^ • (3) 

It is important to notice that (0) is different from the result for on-shell (P 2 = 0) photons 
I9| , due to the different regularization adopted hereQ. This difference will be relevant 
also for the formulation of a model for the parton content of virtual photons, since it is 
part of the perturbatively calculable boundary condition in NLO [pi], [TTH . 

The (Revolutions of the photonic parton distributions are essentially the same for real 
and virtual transverse photons. The inhomogeneous evolution equations are most conve- 
niently treated in the Mellin n moment space, where all convolutions simply factorize, and 



the solutions can be given analytically (see, e.g., [fj3|, |13|). Let us only recall here that the 
distributions f 7 ^ p2 \x, Q 2 ), obtained from solving the inhomogeneous evolution equations, 
can be separated into a 'pointlike' (inhomogeneous) and a 'hadronic' (homogeneous) part 

P^(Q 2 ) = f P P> n (Q 2 ) + fjl^iQ 2 ) . (4) 



In NLO the pointlike singlet solution is schematically given by |2^, [I3fl 



and the usual NLO hadronic solution can be found, e.g, in [ F25| , a, b, d, and U in 
(H) stand for certain combinations of the photon-parton splitting functions and the QCD 
/3-function |23|, [TJ, and L = a s (Q 2 )/a s (Q 2 ). 

Let us finish this technical part by quoting the relevant NLO expression for the struc- 
ture function Fp \x,Q 2 ). It should be pointed out that the treatment and expressions 
for P {p2) (x,Q 2 ) (as on-shell transverse partons obeying the usual Q 2 -evolution equations) 
presented above dictates an identification of the relevant resolved f 1<yP 'X — > X' sub-cross 
sections with that of the real photon according to er^ 7(P )x ~> x ' = 0-f <x -* x \ \ n particular, 
the calculation of \x,Q 2 ) requires the same hadronic Wilson coefficients C 2a and 
C 2 , g as for P 2 = 0, 

= £ 2xe 2 h^+ £ (c 2 , q * q^ + C 2 , 9 * g^) + ^e\C 2 \ + F^P , (6) 

q=u,d,s * 

where F^f ^ represents the charm contribution (see Sec. 4) and * denotes the usual Mellin 
convolution. Note that in the DIS 7 scheme, the NLO direct photon contribution C 2jl in 
(0) is absorbed into the evolution of the photonic quark densities, i.e., C lt2 = 



The difference in the QPM expressions for F^ P ^ between real and virtual photons (as 
pointed out below Eq. (j^)), i.e., in the expressions for C 2tl , is then accounted for by a 
perturbatively calculable boundary condition for q l( - p ' in NLO |13|]. The LO expression 
for Fp is obviously entailed in (|) by dropping all Cj )7 . Finally, it should be noted 



that Fp ^ is kinematically constrained within JTTJ| < x < (1 + P 2 /Q 



}2\-l 



2 Note that F]} F ' is independent of the regularization adopted for calculating "f*(Q 2 )j(P 2 



qq. 
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3 Comparison of different theoretical models 



Let us now briefly highlight the main features of the available theoretical models for the 
parton densities of virtual photons: 

GRS (Gliick, Reya, Stratmann) [|13|| : The GRS distributions provide a straightfor- 



ward and simple extension of the phenomenologically successful GRV photon densities 
2~3| to non-zero P 2 in LO and NLO. As for the GRV densities, the NLO boundary con- 



ditions are formulated in the DIS 7 factorization scheme, originally introduced for P 2 = 
to overcome perturbative instability problems arising in the conventional MS scheme for 
large values of x (see p3| for details). At the low input scale Q = \i ~ 0.5 GeV, uni- 
versal for all 'radiatively generated' GRV distributions (proton, pion, and photon), the 
parton densities of real photons are solely given by a simple VMD-inspired input in LO 
and NL0(DIS 7 ). All one needs to fully specify the distributions for P 2 ^ is a simple, 
physically reasonable prescription which smoothly interpolates between P 2 = (region 
(I)) and P 2 > A 2 QCD (region (III)). This may be fixed by [0 

f^\x, Q 2 = P 2 ) = v(P 2 )fn { o P n -lr t (x, P 2 ) + [1 " V(P 2 )] CPix, P 2 ) (7) 

with P 2 = max(P 2 ,/z 2 ) and f]{P 2 ) = (1 + P 2 /m 2 p )~ 2 where m p refers to some effective 
mass in the vector-meson propagator. Note that the ansatz (|?p implies that the input 
parton distributions are frozen at the input scale fi for real photons for < P 2 < fi 2 such 
that the only P 2 dependence in region (II) stems from the dipole dampening factor i](P 2 ). 

In NL0(DIS 7 ) the perturbatively calculable input /Z^ \x,P 2 ) in Eq. (0) is determined 



by the QPM box result (^); in LO it vanishes (see |13| for details). Since almost nothing 
is known experimentally about the parton structure of vector mesons, the VMD-like non- 
perturbative input is simply taken to be proportional to the GRV pion densities f n |25 



fnl-perti^P 2 ) = * (AlTOl/ f*) X 



T(X,P 2 ) , P 2 >/X 2 

r(x,fi 2 ) , o<p 2 </i 2 



where /z, k, f p are specified in p4 | 



The resulting w-quark and gluon distributions u 7( ^ p2 ^(x,Q 2 ) and g^ p2 \x,Q 2 ), respec- 
tively, are shown in Fig. 1 for Q 2 = 10 GeV 2 and some representative values of P 2 . With 
g~ l<yP ' being the same in the DIS 7 and MS scheme, the shown DIS 7 results for u 1 ^ > can 
be easily transformed to the conventional MS scheme [p3j |24|, [13| . As can be inferred from 
the purely perturbative (rj = 0) contributions, the non-perturbative components, entering 
for rj ^ in Eq. (0), are non-negligible and partly even dominant (especially for x < 0.01). 
It turns out |L3| that only for unexpectedly large P 2 ^> m 2 , say, P 2 (-C Q 2 ) larger than 



about 10 GeV 2 , the perturbative component starts to dominate over the entire x range 
shown (cf., e.g., Fig. 13 in [TJ for P 2 = 20, 100 GeV 2 and Q 2 = 1000 GeV 2 ). 

The precise form of r](P 2 ) in Eq. (|?p clearly represents, apart from fn!m-pert itself, the 
largest uncertainty in this model and has to be tested by future experiments. The only 
measurement of the virtual photon structure in 7*(Q 2 )7(P 2 ) DIS available thus far P?| , is 



7 ( p2\ „ 

compared in Fig. 2 with the NLO GRS prediction for FJ^ = P 2 + |Pl, the combination 
measured effectively by PLUTO [^| (see also Sec. 5). Due to the poor statistics of the 
data and the rather limited x range, the resummed NLO result cannot be distinguished 
from the naive, not resummed QPM result (dotted curve). 




10-1 10-3 10-2 io-l 1 0° 10-4 10-3 10-2 io-l 100 



Figure 1: GRS |T3|] LO and NL0(DIS 7 ) predictions for the w-quark and gluon content of a 
virtual photon for Q 2 = 10 GeV 2 and various fixed values of P 2 (GeV 2 ). For comparison, the 
LO and NLO GRV parton distributions of the real photon (P 2 = 0) []24|] are shown as well. 
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Figure 2: NLO GRS predictions for F]jp = F 2 + \F L |TJ. The data points are taken from 
PLUTO [p6[| . The purely perturbative results correspond to r\ = in Eq. (0). 



Finally, it should be noted that in the GRS approach heavy quarks do not take part 
in the Q 2 -evolution, i.e., there is no 'massless' photonic charm distribution [13]. Heavy 
flavors can only be produced extrinsically, and their contributions have to be calculated 
according to the appropriate massive sub-cross sections (see also Sec. 4). The LO GRS 
distributions are available in parametrized form for P 2 < 10 GeV 2 and Q 2 > 5P 2 p[. 



SaS (Schuler and Sjostrand) ||14| , 15]: The starting point for their analysis are also 



some well-established LO sets of parton densities for real photons |L4j , SaS ID and SaS 2D 



corresponding to two rather different assumptions about the non-perturbative hadronic 
inputf]. The SaS ID set has a similarly low input scale Qo ~ 0.6 GeV as in the GRV [^4]] 
and GRS (13| analyses, but instead of simply relating the VMD input distributions to 
that of a pion, a fit is performed to the coherent sum of the lowest-lying vector meson 
states p, lu, 4>. For the SaS 2D set a 'conventional' high input scale Q = 2 GeV is used at 
the expense of two additional fit parameters, one characterizing the necessary additional 
'hard' component for the quark input at larger values of Q , the other models the effect of 
additional vector meson states beside the ones already taken into account in the SaS ID 
set. The shape of the SaS gluon densities are entirely fixed by theoretical estimates, no 
direct or indirect constraints from direct-photon production data in 7rp collisions as in the 
GRV analysis [^4 have been imposed. Contrary to GRS, heavy flavors are included as 



massless partons in the photon above the threshold Q > m 2 , but when calculating/fitting 
the available F%{x, Q 2 ) data, the massive 'Bethe-Heitler' cross section for 7*7 — > cc is used 
instead, which is, however, not entirely consistent due to double counting. 

The extension to non-zero P 2 is based on the fact that the n moments of the photon 
densities parton can be expressed as a dispersion-integral in the mass k 2 of the 7 — > qq 



fluctuations, which links perturbative and non-perturbative contributions |TJ, [15|. Having 
assumed some model-dependent weight function for the dispersion-integral, and after 
associating the low-k 2 part with some discrete set of vector mesons (as for P 2 = 0), one 
arrives at their final expression for the parton densities of virtual photons 0, 

P ip2) (x, Q 2 ) =J2 ^ h?Jpal Q\ P 2 ) + £ ~e\ r^rP^ix, Q\ k 2 ) , 

y fv l m v + P 1 ^ 7T q Jp2 k 2 

(9) 

where in the perturbative contribution the suppression factor [k 2 /(k 2 + P 2 )] 2 has been 
substituted by an effective lower cut-off P 2 for the integration. Both components f 1,v and 
fi-*qq j n (g) integrate to unit momentum. As in the GRS model above, the VMD part 
contains a dipole suppression factor, which dampens all non-perturbative contributions 
with increasing virtuality P 2 . Contrary the the GRS approach, several different choices for 
the input scale P 2 have been studied apart from P 2 = max(<5o,-P 2 ) The differences 
between all these procedures can be viewed as a measure for the theoretical uncertainty 
within this approach (see, e.g., Figs. 1 and 2 in |T3[). It should be also noted that for some 
more complicated choices for P 2 , the photonic parton densities obey evolution equations 
different from those of the real photon, e.g., the inhomogeneous term can be modified by 
a factor Q 2 j (Q 2 + P 2 ) ||15|| . All different sets of distributions are available in parametrized 
form |L4|, [L5| but, for the time being, the SaS analysis is restricted to LO only. 

Fig. 3 compares the LO GRS with the LO SaS ID and SaS 2D distributions (choosing 
P 2 = max(Qo, P 2 )) for Q 2 = 10 GeV 2 and two P 2 values relevant for future LEP2 mea- 
surements ^7^. As can be seen, the SaS ID results, which refer to a equally low input scale 
as used in GRS, and the GRS densities are rather similar, at least for the w-quark, whereas 
the SaS 2D (quark) distributions are sizeably smaller in this kinematical range, mainly 
due to the higher input scale. For smaller values of x (not shown in Fig. 3) the GRS 
densities rise more strongly than the SaS densities, due the different non-perturbative 

3 The additional SaS 1M and SaS 2M sets jl4| are theoretically inconsistent, as the LO evolved densities 
are combined with the NLO scheme-dependent photon-coefficient function C*2 i7 in the calculation of 
in LO. These sets should not be used in phcnomcnological analyses. 
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Figure 3: Comparison of the LO GRS predictions QT3J] for the w-quark and gluon content of 
virtual photons with the SaS ID and SaS 2D results for P 2 = max(Ql,P 2 ) [jl4[| at Q 2 = 
10 GeV 2 and for two values of P 2 . 

input. However, with increasing P 2 all differences get, of course, more or less washed out, 
since one approaches the purely perturbative domain and the differences in the treatment 
of the non-perturbative component become negligible. 

DG (Drees and Godbole) The aim of this analysis is to estimate the influence 

of the high-P 2 tail in the photon flux in untagged or anti-tagged events, i.e., to study 
the impact of virtual photons in a sample of almost real photons. This effect should be 
taken into account if one tries to extract the parton densities of real photons from such 
measurements. To perform a quantitative analysis, DG provide a simple interpolating 
multiplicative factor r which can be applied to any set of distributions of real photons. 
Several different forms for r are studied in 



and one of the main alternatives is 



1 - 



ln(l + P 2 /P c 



ln(l + Q 2 /P 2 



(10) 



where P c denotes some typical hadronic scale like, for instance, the p mass. The factor r 
is applied to all quark flavors, but the gluon is expected to be further suppressed, since 
it is radiated off the quarks IT! 



q 



^ p2 \x,Q 2 )=rq\x,Q 2 ) , g^ p2 \x,Q 2 ) 



r 2 g\x,Q 2 ) 



;n) 



For typical 77 experiments with Q 2 ~ 10 GeV 2 in an untagged situation, virtual photon 
effects then suppress the effective photonic quark and gluon content by about 10 and 15%, 
respectively |16| . 

Obviously, the above ansatz (O) does not change the x shape of the distributions, 



which would require more complicated forms for r [[U|. The recipe (|TT[) does not appear 
to be well suited for QCD tests of the virtual photon content: on the one hand the densities 
in (|TT|) are not a solution of the inhomogeneous evolution equations, and on the other 
hand there is no dipole power suppression factor for the non-perturbative VMD part of 
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the virtual photon densities as in GRS [13 



[15 1, i.e., the approximation (|TTD 



or SaS 

can only be applied in the large x region where the perturbative pointlike part dominates, 
whereas at smaller x it may grossly overestimates the densities for increasing virtuality 
P 2 . A similar strategy as in Eqs. (|ToD and ([TT| ) has been used in [18] including, however, 
a power suppressed VMD part. 



4 Treatment of heavy flavors 

The question of how to treat heavy flavor (m q 3> Aqcd) contributions to structure func- 
tions and cross sections in the most appropriate and reliable way has attracted a consid- 



erable amount of interest in the past few years |2T], 0, [28]. This was mainly triggered 
by the observation that the charm contribution to the DIS proton structure function Ff 
amounts to about 20 — 25% in the small x region covered by HERA. In case of the photon 
structure function Fj in (||) , effects due to charm are sizeable also in the large x region 
due to the existence of the direct/pointlike component, such that a proper treatment is 
even more important here. 

There are two extreme ways to handle heavy flavors: one can simply include heavy 
flavors as massless partons in the evolution above some threshold Q 2 > m 2 , or one can 
stick to a picture with only light partons in the proton/photon. In the latter case, heavy 
flavors do not participate in the evolution equations at all and can be produced only 
extrinsically. In case of F^f in (Jfty, two different contributions have to be taken into 
account. Firstly, the direct 'Bethe-Heitler' process Y(Q 2 )l(P 2 ) ~^ and secondly the 
resolved contribution 7*(Q 2 )(7 7 ^ P ' — > cc. For real photons these cross sections are known 



up to NLO and can be found in p9|. In |30[ it was shown that the two contributions 



are separated in the variable x: for large > 0.05, the direct process dominates, 

whereas for x < 0.01 the dominant contribution stems from the resolved part. 

The fully massless treatment has been abandoned recently in all existing modern sets 
of proton densities, simply because it does not exhibits the correct x and Q 2 dependent 
threshold behavior. A potential problem with the massive treatment are possibly large 
logarithms in the relevant sub-cross sections far above the threshold, which might be 
indicative for resummations, i.e., for introducing a 'massless' heavy quark distribution. 
Therefore various 'unified' prescription were proposed recently f22| which reduce to the 



massive results close to, and to the massless picture far above threshold. For the time 
being these studies have been performed only for the proton densities and not in the 
context of photons. 

However, as mentioned in Sec. 3, GRS [13| prefer to stick to the fully massive frame- 



work, similarly to the case of the GRV proton densities [[31]]. In each case this is motivated 
by the observation that all relevant fully massive production mechanisms appear to be 
perturbatively stable, at least for all experimentally relevant values of x and Q 2 (see, e.g., 
Refs. pl|, p0|). Moreover, all theoretical uncertainties, in particular the dependence on 
the factorization scale appear to be well under control, leading to the conclusion that 
there is no real need for any resummation procedure. It should be mentioned that the 
relevant expressions for 7*(Q 2 )7(-P 2 ) — > cc for non-zero P 2 are available only in LO so far 
|32|| , hence a study of the perturbative stability cannot be performed here yet. 

Clearly, the fully massive treatment is much more cumbersome and inconvenient than 
the massless or 'unified' framework when calculating, for instance, jet production cross 
sections in ep or 77 collisions. One cannot simply increase the number of active flavors by 







one unit and use a c 7< - p2 )(x, Q 2 ) distribution. Instead one has to calculate the relevant sub- 
cross sections with massive quarks for the final state configuration under consideration, 
which is much more involved and time-consuming in numerical analyses. However, for 
large-pr jet production in LO (m q /pT <C 1), for instance, one can simply approximate 
the relevant massive cross sections by their massless counterparts, neglecting of course all 
massless contributions with a 'heavy' quark in the initial state. 



5 Measuring 7*-PDF's in e + e and ep reactions 

Since there are several dedicated contributions which discuss recent experimental progress 
or future prospects ||, [F], |33| , we can be fairly brief here and concentrate only on topics 
not covered elsewhere. 

Let us first of all delineate the x, Q 2 , and P 2 ranges covered by LEP2 and the expected 
statistical accuracy for the virtual photon structure functions measurements [p7[| . Because 



of its higher energy and integrated luminosity, LEP2 can provide improved information 



from double tagged events as compared to the not very precise results from PLUTO |2C 
shown in Fig. 2. The most important double tag sample at LEP2 is expected to come 
from events with Q 2 > 3 GeV 2 and 0.1 < P 2 < 1 GeV 2 . For typically expected 500 pb~ x 
of data collected, about 800 of such events will be seen, covering 3 • 10 -4 < x < 1 and 
3 < Q 2 < 1000 GeV 2 (Z7|. The yield of events with both Q 2 and P 2 > 3 GeV 2 seems to 



be too small for a meaningful analysis. 

Fig. 4 shows the virtual photon structure function F^ P \x, Q 2 ) in LO as predicted by 
the SaS ID, SaS 2D (using P 2 = max(Q 2 , P 2 )), and GRS models in two bins for P 2 and Q 2 . 
The error bars indicate the statistical precision expected for each x bin using the SaS ID 
densities (similar for the SaS 2D and GRS distributions). A measurement of F^ P (x, Q 2 ) 
within these bins, as distinct from the real (P 2 = 0) photon structure function (illustrated 
by the solid curves for SaS ID) should be possible at LEP2 and could be compared to the 
different model predictions, which turn out to be rather similar in the accessible x, Q 2 , 
and P 2 bins, except for the smallest x bins and P 2 — > ((P 2 ) = 0.2 GeV 2 ). The latter 
differences are of course related to the present ignorance of the P 2 = distributions for 



x — > 0, i.e., whether they either steeply rise as in case of GRV [23] or show a rather flat 
x — > behaviour as, e.g., in case of SaS 2D |l4j, [15 . 



However, apart from the experimental challenge there is an additional complica- 
tion already noticed by PLUTO [^6|: what is directly measured is, of course, not 

F] but the / y*(Q 2 ) / y(P 2 ) DIS cross section, which can be schematically expanded as 
0" 7 .(q2) 7 ( P 2) = o"tt + £iO~lt + £20~tl + where L and T denote longitudinal and 

transverse polarization, respectively, of the probe and the target photons, and are the 
L/T 7-flux ratios. For PLUTO e\ ~ e 2 ^ 1 y <S 1) and assuming that ctll — 
and o~lt = o~tl, as for the QPM expressions for Jt l(Q 2 )1t,l(P 2 ) —> qq for vanishing con- 
stituent quark masses []3~2]| , one arrives at the combination |26| cr 7 . 7 ~ F 2 + 3/2F L = F^ P 



eff 



effectively measured by PLUTO (cf. Fig. 2). Hence, strictly speaking such measurements 
cannot be directly related to the densities f 1<yP \x,Q 2 ), since only transverse (T) virtual 
photons are described by the GRS, SaS, and DG models. Furthermore, in the QPM model 
|32|1 it turns out that the contribution due to longitudinal target photons is rather size- 



able at large x even for P 2 /Q 2 -C 1, contrary to the expectation that transverse photons 
should dominate in this region. Clearly, more work is required here for a meaningful in- 
terpretation of any future results from LEP2, possibly including also studies of the parton 
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Figure 4: Expectations for the statistical accuracy of the virtual photon structure measurement 
at LEP2 in two different P 2 and Q 2 bins using the SaS ID distributions p| . The SaS ID 
predictions for P 2 = and the results for the GRS [JT3]] and SaS 2D models are shown as lines 
for comparison. The upper (lower) curves for GRS and SaS 2D refer to P 2 = 0.2 (0.5) GeV 2 , 
respectively. The figure is taken from p7| . 

content of longitudinal photons which have not been carried out so far. 

In ep collisions (di-)jet production is certainly the best tool to decipher the parton 
structure of virtual photons, and a lot of experimental and theoretical progress was re- 
ported at the workshop j|, [/]]. There should be a hierarchy between the hard scale ji 2 



Q 2 in e + e ) at which the virtual photon is probed (typically fij 



0(P 2 + Vhet) 



in case of jet production) and the photon virtuality P 2 (see footnote 1). Exactly in 
this kinematical domain an excess in the dijet rate was observed by HI [[|, which can be 
nicely attributed to a resolved virtual photon contribution, in accordance with all existing 
models for the •* described in Sec. 3. 

Recently, similar studies were extended to the production rate of forward jets |3"B| , 
which is regarded as a test in favor of the BFKL dynamics |34|. Indeed the usual DGLAP 
(direct photon induced) cross section falls short of the data by roughly a factor of two |35| . 
In it was demonstrated, however, that the inclusion of the resolved virtual photon 
component removes this discrepancy, and the full DGLAP results are then in even better 
agreement with data than the BFKL results, in particular for the two to one forward 
jet ratio |]35|| . However one should be cautious to jump to any conclusions. In order to 
suppress the phase space for the DGLAP evolution, the p\ of the forward jet is required 
to be of the same size as the virtuality P 2 of the photon, hence there is no real hierarchy 
between the hard scale ji 2 = 0{P 2 +Pr,jet) an d P 2 an d thus no large logarithm ln/z 2 /P 2 



which is resummed in / 



Naively one would therefore expect only a small resolved 



contribution, as was also observed in the HI jet analysis [Q or in the theoretical studies |7[] 
for P 2 — > fi 2 , rather than a gain by about a factor of two. Hence the kinematics of forward 
jets seems to be very subtle (for instance, the virtual photon content is only probed at 
large momentum fractions x 7 ), and presumably the theoretical uncertainties due to scale 
variations and changes in the model for the f^ p •* (in particular, of the input scale P 2 ) 



1 1 



are of the same size as the resolved photon contribution itself. More detailed studies are 
clearly required here. Furthermore, it should be kept in mind that all BFKL results so 
far are based only LO parton level calculations. It seems, however, that the forward jet 



kinematics is not suited to distinguish between BFKL and DGLAP at HERA [ 33 1 
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